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Recently developed strong-coupling theory open up the 
possibility of treating quantum-mechanical systems with 
hard-wall potentials via perturbation theory. To test the 
power of this theory we study here the exactly solvable quan- 
tum mechanics of a point particle in a one-dimensional box. 
Introducing an auxiliary harmonic mass term m, the ground- 
state energy E^ ' can be expanded perturbatively in powers 
of 1/md, where d is the box size. The removal of the infrared 
cutoff 771 requires the resummation of the series at infinitely 
strong coupling. We show that strong-coupling theory yields a 
fast-convergent sequence of approximations to the well-known 
quantum-mechanical energy E^°' = 7r 2 /2d 2 . 

1. Variational perturbation theory H permits us to 
convert divergent weak-coupling expansions into conver- 
gent strong-coupling expansions. In particular, a con- 
stant strong-coupling limit of a function can be evalu- 
ated from its weak-coupling expansion with any desired 
accuracy. As an important application, this has led to a 
novel way of calculating critical exponents without using 
the renormalization group |g| . 
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Given this theory, new classes of physical systems be- 
come accessible to perturbation theory. For instance, the 
important problem a the pressure exerted by a stack of 
membranes upon the enclosing walls [p[ has now become 
calculable analytically with the help of perturbation the- 
ory. For a single membrane, this has already be done 
successfully [f|||. Realistic physical problems have usu- 
ally the disadvantage that the maximally accessible or- 
der of perturbation theory is quite limited. If we want 
to gain a better understanding of the convergence of the 
successive approximations as the order goes to infinity it 
is useful to study a system where the result is known 
exactly. This will be done in the present note for a 
quantum-mechanical point particle in a one-dimensional 
box. The ground state energy of this system is known 
exactly, E^ — n 2 /2d 2 (in natural units), where d is size 
the of the box. We shall demonstrate how this result 
is found via strong-coupling theory from a perturbation 
expansion, thus illustrating the reliability of the earlier 
membrane calculations 

2. The partition function of a particle in a box is given 
by the euclidean path integral (always in natural units) 



= (Vu(t)e^ dt ^\ 



(1) 



where the particle coordinate u(t) is restricted to the 
interval —d/2 < u(t) < d/2. Since such a hard-wall 
restriction is difficult to treat analytically in the path 
integral (||) , we make the hard- walls soft by adding to the 
euclidean action E in the exponent of (0) a potential term 
diverging near the walls. Thus we consider the auxiliary 
euclidean action 

E= l -Jdt{[du{t)f + V{u{t))}, (2) 

where V(u) is given by 
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y(u) = TU tan d 



On the right-hand side we have introduced a parameter 
g^n 2 /d 2 . 

3. The expansion of the potential in powers of g can 
now be treated perturbatively, leading to an expansion 
of Z around the harmonic part of the partition function, 
in which the integrations over u{t) run over the entire 
w-axis and yield 
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For L — > oo, the exponent gives a free energy density 
/ = — L _1 log Z equal to the ground state energy of the 
harmonic oscillator 
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The treatment of the interaction terms can be organized 
in powers of g, and give rise to an expansion of the free 
energy with the generic form 
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The calculation of the coefficients a,k in this expansion 
proceeds as follows. First we expand the potential in (g) 
to identify the power scries for the interaction energy 
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The interaction terms J dt [u 2 (t)] k+1 and their products 
are expanded according to Wick's rule into sums of prod- 
ucts of Wick contractions representing harmonic two- 
point correlation functions 

(u(ti)u(ta)) = / 7ZZ -TT—7T = K~. ■ ( 9 ) 
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Associated local expectation values are (u 2 ) = l/2ui, and 
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where the last integral is calculated using dimensional 
rcgularization in which J dkk a = for all a. The Wick 



contractions are organized with the help of Feynman di- 
agrams. Only the connected diagrams contribute to the 
free energy density. The graphical expansion of free en- 
ergy up to four loops is 
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Note different numbers of loops contribute to the terms 
of order g n . The calculation of the diagrams in Eq. (|l] 
is simplified by the factorization property: If a diagram 
consists of two subdiagrams touching each other at a 
single vertex, the associated Feynman integral factorizes 
into those of the subdiagrams. In each diagram, the last 
t-integral yields an overall factor L, due to translational 
invariance along the i-axis, the others produce a factor 
1/oj. Using the explicit expression ( |lO| ) for the lines in the 
diagrams, we find the following velues for the Feynman 
integrals: 
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Adding all contributions in (|11|), we obtain up to the 
order g : 
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which has the generic form ^). We can go to higher 
orders by extending the Bender- Wu recursion relation 
for the ground-state energy of the quartic anharmonic 
oscillator as follows: 



2jC„ J = (j + l)(2j + l)C„ J - 
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(14) 



After solving these recursion relations, the coefficients a,k 
in (H) are given by au = (— l) k+1 Ck,i- For brevity, we 
list here the first sixteen expansion coefficients for /, cal- 
culated with the help of the algebra program REDUCE: 
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4. We are now ready to calculate successive strong- 
coupling approximations to the function g(g) and suty 
the convergence behavior as the order grows large. Ac- 
cording to the general theory in [g]jl]] , the Nth order ap- 
proximation to the strong-coupling limit /* of the series 
f(g) is found by replacing, in the truncated series after 
the Ath term /jv(<7), the frequency to by the identical 
expression \f$l 2 + gr 1 where r = (oj 2 — Q 2 )/ g which is, 



however, treated for a moment as an independent vari- 
able, whereas fHs a dummy parameter. Then the square 
root is expanded binomially in powers of g, and /at(<7) 
is reexpanded up to order g . After that, r is replaced 
by its proper value. In this way we obtain a function 
/jv(3) f^) which depends on O, which thus becomes a vari- 
ational parameter. The best approximation is obtained 
by extremizing /jv(<7,f2) with respect to w. The strong- 
coupling limit is obtained by taking g to infinity which 
is equivalent to setting u> = 0. In this limit, the optimal 
n will grow proportionally with g, so that g/Ct = 1/c is 
finite, and the variational expression /jv(g, fi) becomes 
a function of /at(c). For oj = 0, the above reexpan- 
sion amounts simply to replacing each power uj n in each 
expansion terms of /at (g) by the binomial expansion of 
(1 — 1)~™/ 2 truncated after the (TV — n)th term. The first 
nine variational functions /at(c) are listed in Table 1. 

TABLE I. First eight variational functions /jv(c). 
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The functions /at (c) are minimized starting from f% (c) 
and searching the minimum of each successive /3(c), 
/3(c),... nearest to the previous one. The functions 
/at(c) together with their minima are plotted in Fig. fy. 




FIG. 1. Variational functions /jv(c) up to iV = 16 are 
shown together with their minima whose y-coordinates ap- 
proach rapidly the correct limiting value 1/2. 

The minima lie at 

{N, f$ iu ) = (2,0.466506), (3,0.492061), (4,0.497701), 
(5,0.499253), (6,0.499738), (7,0.499903), 
(8, 0.499963), (9, 0.499985), (10, 0.499994), 
(11,0.499998), (12,0.499999), (13,0.5000), 
(14,0.50000), (15,0.50000), (16,0.5000). (16) 

They converge exponentially fast against the known re- 
sult 1/2, as shown in Fig. 0. 
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FIG. 2. 



5. The alert reader will have noted that the expansion 



coefficients ( |15[ ) possesses two special properties: First, 
they lack the factorial growth at large orders which would 
be found for a single power [u 2 (£)] fc+1 of the interaction 
potential . The factorial growth is canceled by the spe- 
cific combination of the different powers in the interaction 
(M) , making the series (ra) convergent inside a certain cir- 
cle. Still, since this circle has a finite radius (the ratio 
test shows that it is unity), this convergent series cannot 
be evaluated in the limit of large g which we want to do, 
so that variational strong-coupling theory is not super- 
fluous. However, there is a second remarkable property 
of the coefficients (J15J) : They contain an infinite num- 
ber of zeros in the sequence of coefficients for each odd 
number, except for the first one. We may take advan- 
tage of this property by separating off the irregular term 
a±g = g/A = 7r 2 /4eP, setting a — g 2 /Aui 2 , and rewriting 
/ as 

N 
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Inserting the numbers (|15|), the expansion of h(a) reads 
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Here we can go directly to the strong-coupling limit a — > 
oo corresponding to d — > 0, where we recover the exact 
ground-state energy E^ * 1 — it 2 /2d 2 . 

5. The energy ( ^9[ ) can of course be obtained directly 
by solving the Schrodinger equation associated with the 
potential (R), 
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where we have replaced u — » dx/ir and set w 2 d 4 /7r 4 = 
A(A — 1), so that 
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Equation ( p0| ) is of the Poschl- Teller type and has the 
ground state wave function pf] 



V>o (x) = const x cos A x , 



(22) 



with the eigenvalue tt 2 E^ jd 2 = (A 2 — 1)/2, which agrees 
of xourse with Eq. (In 

If we were to apply the variational procedure to the 
series h(a)/y/a in / of Eq. (JlJ), by replacing the factor 



(18) 



We now realize that this is the bibomial power series 
expansion of \/\ + a. Substituting this into ( |17J ) , we find 
the exact ground state energy for the euclidean action 



\j{jj 2n contained in each power a n by fi — \/£l 2 + ra 
and reexpanding now in powers of a rather than g, we 
would find that all approximation /ijv(c) would posses a 
minimum with unit value, such that the corresponding 
extremal functions /n (c) yield the correct final energy in 
each order TV. 
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